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AN ANALYSIS OF LATERAL STABILITY IN POWER-OFF FLIGHT 
WITH CHARTS FOR USE IN DESIGN 

By Charles H. Zimmerman 



SUMMARY 

The aerodynamic and mass factors governing lateral 
stability are discussed and formulas are given for their 
estimation. Relatively simple relationships between the 
governing factors and the resulting, stability characteristics 
are presented. A series of charts is included with which 
approximate stability characteristics may be rapidly . 
estimated. 

The effects of the various governing factors upon the 
stability characteristics are discussed in detail. It is 
pointed out that much additional research is necessary both 
to correlate stability characteristics with riding, flying, and 
handling qualities and. to provide suitable data for accurate 
estimates of those characteristics of an airplane while it is 
in the design stage. 

INTRODUCTION 

The lateral stability of airplanes has been the subject 
of considerable mathematical treatment and many- 
theoretical analyses. (See references.) The main as- 
pects of the problem are therefore well known to stu- 
dents of the subject. Use of the mathematical theory 
in design is, however, limited by practical difficulties in 
its application. Determination of numerical values for 
certain of the aerodynamic quantities is difficult and 
the results are uncertain. The required calculations 
are extensive and must be carefully made to avoid erro- 
neous and confusing results. 

In this report lateral stability will be discussed and 
analyzed in a way that, it is believed, will aid in the 
acquisition of a working knowledge of the subject with- 
out long and intensive study. The classical equations 
have been simplified as much as seems consistent with 
reasonable accuracy to permit rapid estimation of the 
stability characteristics. Also included is a series of 
charts designed to facilitate the rapid estimation of the 
approximate lateral-stability characteristics of airplanes 
throughout the normal-flight range. It is hoped that 
these charts, together with those on longitudinal sta- 
bility presented in reference 1, will aid in putting the 
estimation of the complete stability characteristics on 
a practical basis. 

The material is presented in the following order: (1) 
A discussion of the aerodynamic and mass factors that 



govern the uncontrolled motion of the airplane together 
with formulas for estimating these factors; (2) formulas 
for estimating the stability characteristics of the uncon- 
trolled motion having given the governing factors; (3) 
charts for the rapid estimation of stability character- 
istics; (4) a discussion of the effects of the governing 
factors upon the stability characteristics; (5) comments 
and suggestions for future study; (6) a brief derivation 
of the classica| stability formulas (appendix I); (7) an 
accurate semigraphical method for solving biquadratics 
with a useful approximation based on this method 
(appendix II); and (8) a list of symbols and their 
definitions (appendix III). 

FACTORS GOVERNING STABILITY 

. Both theory and experiment indicate that, with 
certain exceptions, the uncontrolled motion of an air- 
plane can be divided into two independent phases. 
One 'phase includes components of the motion that do 
not displace the plane of symmetry of the airplane 
from the plane with which it coincides during the steady 
motion. Stability of this part of the motion is termed 
"longitudinal stability." The other phase of the com- 
plete motion includes all components that do displace 
the plane of symmetry. This phase of the motion is 
called "lateral motion" and its stability characteristics, 
"lateral stability." Although, in the past, reference 
has frequently been made to directional stability as 
distinguished from rolling stability (also called "lateral" 
stability), both theory and experiment indicate that 
no such division is physically possible for the conven- 
tional airplane. 

The uncontrolled motion of an airplane quite ob- 
viously depends upon the aerodynamic forces and 
moments arising from any deviation from a steady 
state together with the inertial forces and moments 
accompanying the accelerations coupled with the 
deviations. The lateral motion is zero in steady 
flight on a straight course. The components of lateral 
motion in unsteady flight are a linear velocity v along 
the Y axis (see appendix I and fig. 1) and angular 
velocities p and r about the X and Z axes, respectively. 
The forces and moments governing lateral motion 
therefore arise from the aerodynamic reactions to the 
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velocities v, p } and r (in the theoretical treatment aero- 
dynamic reactions are assumed to be unaffected by 
accelerations) and the inertial reactions to the accelera- 
tions dvfdt,g sin <f> cos 7, g sin ^ sin 7, dp/«#,.and dr(dt, 
where <j> is the angle of roll, 7 is the angle of the flight 
path, and ^ is the angle of yaw. 

For convenience the components of the reactions 
referred to the coordinate axes are used rather than 
the resultant reaction. It appears, then, that a velocity 
v should result in a side force AY, a rolling moment AL, 
and a yawing moment AN. Similarly there will be 
A's of Y, L, and N corresponding to the rolling and 
yawing velocities p and r. The basis for the classical 
theory of stability is that the algebraic sum of the 
values of AY (for example) for a unit value of v when p 
and r are zero, for a unit value of p when v and r are 
zero, and for a unit value of r when v and p are zero is 
equal to the value of AY when the total motion is the 
resultant of coexisting unit values of v, p, and r. It is 
further assumed that a reaction AY due to a disturb- 




Figube 1.— Angular and vertical relationships in flight, power ofi. 

ance of velocity v is directly proportional to the magni- 
tude of v, that ia AY—v^- This assumption is ad- 
mittedly an approximation but is valid, in general, for 
small values of the velocities of the disturbance. On 
this basis the aerodynamic reaction AY to a lateral 
disturbance is 

dv'^dp^^dr 
and similar expressions exist for AL and AN. 

As a matter of convenience it has been found desir- 
able to express the derivatives dYjdv, dL/dv, etc., in 
terms of the nondimensional coefficients C Y , <7„ and C n 
where 



N 



In order to make the treatment entirely nondimensional, 
it is convenient to consider the ratios v/V, pb/2V, and 
rb/2V rather than v, p, and r. For small values v/V 
is equal to 0, where 0 is the angle of sideslip (in radians), 
and pb/2V is the difference (in radians) between the 
angle of attack at the center of gravity and the angle 
of attack at the wing tip. Since the velocity at the 
wing tip is V+rb/2 the value rb/2V is the ratio of the 
portion of the velocity at the tip due to rotation to the 
velocity at the center of gravity. Expressed in this 
way, the lateral-force coefficient due to lateral motion is 



dOy , pb dCy 



2Vd 



pb 



rbdCy 



'2V 



and similar expressions exist for Ad and AC n . 

Since dC Y /d£p and dC Y /d^a,re small, they are gen- 
erally neglected, leaving the following aerodynamic 
factors to be considered: 

1. Those depending on sideslip: dC y /dp, dCJdp, and 
dC n /d0. 

2. Those depending on rolling velocity: dCJd^ and 
dCJd^- 

3. Those depending on yawing velocity: ddjd^. and 



dOJd, 



rb 



2V 



In addition to the aerodynamic factors, others that 
depend on the amount and the distribution of the mass 
of the airplane must be considered. The important 
mass factors, expressed nondimensionally, are /x, bfk x , 
and b}k z . The relative density factor p is equal to 
mjpSb and may be considered as being proportional to 
the ratio of the mass of the airplane to the mass of air 
influenced by it in traveling one chord length. Under 

standard conditions ff= 13 - 1 W'S) . 

0 

AERODYNAMIC FACTORS 

lateral force due to sideslip— The rate of change of 
lateral-force coefficient with angle of sideslip dC Y /dp can 
be accurately determined only by measurement in a 
wind tunnel. Assuming the wind-tunnel data to have 
been obtained in terms of angle of yaw ^ in degrees, the 
value of dCr/dfi is —57.3 (dC Y /df), since j8 is in radians 
and opposite in sign to In wind-tunnel practice, 
cross-wind force rather than lateral force is usually 
measured. In such cases dC Y /dp can be determined 
from the relationship 



dCr__dC 6 ' 
dfi~dfi^ tD 



(1) 
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(which follows from the fact that C Y =C e cos 0 
— C^sin/S). 

Diehl gives (reference 2, pp. 254-255) an approxi- 
mate, empirical value of 



^— 0 12* 



(2) 



where Z, is the over-all length. This formula is useful 
when wind-tunnel data are not available. 

Rolling moment due to sideslip. — The rate of change 
of rolling-moment coefficient with sideslip dCJdp must 
also be measured in a wind tunnel if accurate values are 
desired. Some systematic research has shown the effect 
of dihedral and tip shape on the value of dCJdp for the 
wing alone (reference 3) but very little is known about 
the effect of fuselage interference. In certain experi- 
ments (data unpublished) a model having a wing with 
no dihedral mounted in a high-wing position gave a 
value of dCt/dp corresponding to 5° of positive dihedral 
for the wing alone. The same model with the wing 
mounted in a low-wing position gave a very erratic 
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Figure 2.— F.?oct of tip shape on rate of change of rolling-moment coefficient with 

sideslip. 



curve of C t against p. The slope of this curve indicated 
zero dihedral effect at zero sideslip. The average di- 
hedral effect up to 30° sideslip corresponded, however, 
to 4° negative dihedral. These tests were in the nature 
of preliminary tests and are unconfirmed but give ample 
evidence of the need for similar additional research. 



In the absence of wind-tunnel tests the value of 
dCJdfi for the wing alone may be computed from the 
relationship 

' f=(iL +r (-°- 012 > m 

where (dC t ldp) T . 0 is the value of dOJdp for the wing 
without dihedral (see fig. 2) and T is the dihedral angle 
in degrees. This formula was developed from data ob- 
tained with wings of aspect ratio 6 and with no taper or 
sweepback (reference 3). Tapering the wing decreases 
the effective dihedral but the decrease is somewhat less 
than would be expected from the geometric proportions 
because of the tendency of the wing lift to be evenly 
distributed along the span. Sweepback is equivalent 
to an increase in dihedral, particularly at high angles 
of attack, but the effect is negligible for small amounts 
of sweepback such as are used in conventional airplanes. 

The wing, including interference effects, is the chief 
source of rolling moment due to sideslip and other parts 
of the airplane can normally be neglected. Vertical- 
fin area displaced from the longitudinal axis contributes 
to dC x \dp but the effect is usually small. If, in a par- 
ticular case, the effect upon the value of dCJdp is de- 
sired for parts having considerable projected side area, 
it can be computed from the relationship 
dCi_S s z R dC Lp 

dp~Sb dp W 
where S P is the projected side area. 
S t the wing area. 

z v , the z coordinate of the center of pressure of 
projected side area. 
* C Lpr the absolute coefficient of force on the pro- 
jected side area. 
In this equation dC Lp jdp must be estimated, taking 
into account the shape of the part and the probable 
interference effects. 

Yawing moment due to sideslip. — The change of 
yawing-moment coefficient with angle of sideslip 
dCJdp depends principally upon the fuselage and the 
vertical-tail area. The contributions of the landing 
gear, interference effects, etc., are small and can gen- 
erally be neglected. The contribution of the wings 
is also small and can be neglected at high or cruising 
speeds but becomes of increasing importance at slower 
speeds (reference 3); the effect due to the wing is an 
increase in dGjdp. The center of pressure upon the 
fuselage is normally well ahead of the center of gravity 
so that the moment due to sideslip is such as to increase 
the sideslip. The magnitude of this unstabilizing 
tendency varies with the' size and shape of the fuselage 
but, on the average, is equal to about one- third the 
stabilizing effect of the vertical tail surfaces. 

For accurate stability calculations it is necessary that 
dOJdp be obtained from wind-tunnel tests at several 
angles of attack by the use of the relation 
dC„/dp= -57. 3 dCJdf. The value of dCJdp can be 
calculated approximately from the relation 
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dCn_ XSt/dOj 



■Kf 



Sb 



(5) 



where 77 , is the tail efficiency. 

I, the distance from the center of gravity to 

the rudder hinge. 
k, the over-all length of the fuselage.. 
S t , the area of the vertical tail surfaces. 
<7ii, the absolute coefficient of force on the verti- 
cal tail surfaces. 
Kp, an empirical constant (reference 2, p. 203). 
S t , the projected side area of the fuselage. 
When using equation (5), it is necessary to estimate or 
assume values of rj t} dCu/dfi, and Kp. For modern 
types of airplane r\ t is about 0.80. The slope of the 
tail-force curve dCn\d$ depends on the aspect ratio of 
the tail and to a certain extent upon the end-plate 
effect of the fuselage and the horizontal surfaces. For 
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Figure 3.— Empirical factor for computing effect of fuselage on rate of change of 
yawing-moment coefficient with sideslip (from flg. 98, reference 2). 



conventional arrangements dCL t jd$ = 2.2 is a good 
average value. Values of Kp as determined by Diehl 
are given in figure 3. From this figure the value of K$ 
can be directly determined from the ratio of the distance 
of the center of gravity back of the nose Xi to the fuselage 
length lx and the ratio of the maximum fuselage depth d 
to the fuselage length. In a number of computations 
made to check the accuracy of formula (5) it was found 
that the results were generally conservative, i. e., the 
estimated value of dCJdfi was smaller than the meas- 
ured value. . The difference arose in most cases from 
the fact that the measured effect of the fuselage was 
smaller than the estimated effect. The measured 
effect of the fuselage apparently varies between zero 

Cf 7 

and the effect calculated as depending upon the 

details of nose shape and fuselage foim. 



Rolling moment due to rolling, — The rate of change 

T)b 

of rolling-moment coefficient with rate of rolling dCi/d^y 

arises from the change of angle of attack along the 
wing. The increment in angle of attack at any span- 
wise distance y from the center of gravity is py/V (in 
radian measure), the increment at the tip being pb/2V. 
If a uniform spanwise distribution of lift and drag be 
assumed, simple integration gives 

dGi^ -idajda + Gp.) 

d2 b 6 

a 2V 

where C Dw is the drag coefficient of the wing alone. 
If an elliptical distribution is assumed, integration gives 

dCi _~(dC L /da+C Dw ) 



d^~ 
a 2V 



8 



Actually, measured values of dCifdj^y&re considerably 

smaller than either formula indicates because of the 
tendency of the lift to equalize itself along the span 
during the rotation. 
In the absence of data obtained from some such device 

as a rolling balance, dOi/d^p- can betaken as —0.40 

for wing arrangements such as are likely to be used on 
conventional airplanes. A survey of test results 
reveals values from —0.35 to —0.47 for plain wings 
and values as high as —0.50 for wings with tip slots. 
It would be expected that rounding the tips or tapering 

the wings would reduce d/7 t / d^p^and such was found to 

be the case for the tests reported in reference 3. On 
the other hand, there is sufficient conflicting evidence to 

vb 

indicate that an attempt to calculate dGi/d^y taking 

into account tip shape and taper, is likely to give a 
result no nearer the true value than is the assumed 
average, —0.40. 

Yawing moment due to rolling. — The rate of change 
of yawing-moment coefficient with rate of rolling 

dC n jd^ arises from the same causes as does ddjd-^- 

Simple integration gives for a rectangular-wing force 
distribution 



d 2% 



and for an elliptical-wing force distribution gives 



dCn_ 

X 2V 
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It will be noticed that the sign of the resulting value of 
t)b 

dCnfd^y indicates that the wing being depressed by the 

rolling motion is accelerated forward by the resulting 
yawing moment. The mistake has frequently been 
made (see reference 4) of assuming that the increase in 
drag of the wing being depressed would result in a 
yawing moment retarding that wing, that is, in a 

positive value of dC n /d^ This reasoning fails to 

take into account the forward inclination with increase 
in angle of attack of the resultant-force vector relative 
to an axis fixed in the wing. 
Wind-tunnel data cannot ordinarily be obtained for 

dC n /d^y because there are but few existing balances 

capable of measuring the yawing moment on a rolling 
model. It is therefore necessary to rely on estimated 

values of dC n /d£p> The empirical relationship 



dG n 

d&- 
a 2V 



mm 



(6) 



has been found (reference 5) to give good agreement 
with measured values below the stall, O h and dC D Jda 
having been obtained from force tests of the wing 
alone; but there is need for further experimental data 
on this factor. 

Rolling moment due to yawing. — The rate of change 
of rolling-moment coefficient with rate of yawing 
rb 

dCi/d^y results from the difference in velocity between 

the wing tips, one wing tip having the velocity V+rbf2 
and the other having the velocity V—rb/2. Simple 

integration gives the value of dd/d^p. as C L {Z assuming 

a rectangular distribution of lift or as Oj4c assuming an 
elliptical distribution. The rolling moment due to 
yawing is of positive sign since a positive rate of yawing 
gives a positive rolling moment. 

It will ordinarily be impossible to measure dC x \d£t 7 

for a particular design because of lack of equipment. 
Either special apparatus for oscillating the model or a 
whirling arm equipped to measure rolling moments is 
required. In the absence of experimental data, the 
computed value must be used. Glauert states (refer- 
ence 6) that 

.rb 4 W 
a 2V 

gives nearly correct values for a rectangular wing. 
Experimental results for the Bristol Fighter, a biplane 
with substantially rectangular wings, however, gave 



dCijd^y as nearly 0^/3 (reference 7). Measured values 

from tests of a biplane model reported in reference 8 
were approximately equal to C L /4 for three wing combi- 
nations. It appears that the assumption that equation 
(7) gives reasonable values is justfied for wings with 
faired or elliptical tips and slight to moderate taper. 
Yawing moment due to yawing— The rate of change 

of yawing-moment coefficient due to yawing dC n jd^ 

results from the change of velocity along the wing and 
the change of sideslip velocity along the fuselage and 
at the tail due to the yawing. On the basis of simple 

integration the portion of dOJd— due to the wing is 

— 0 D JZ for a rectangular distribution and —C D j± for 
an elliptical distribution. In an extension of work by 

Wieselsberger, Glauert shows (reference 6) that dG„/d~ 

is equal to -(0.33 £7^+0.043 C Dl ) for a rectangular 
wing of normal aspect ratio and equal to —(0.25 0 D + 
0.33 C Di ) for an elliptical wing, where C Do and C Di are 
the profile and induced drags, respectively, for the 
wing alone. 

The change in angle of sideslip at the tail due to a 
yawing velocity r is rl/V. The theoretical value for the 
vertical tail is 

dO n ■ P8 t dC Lt 
jrb ~ ZVt b 2 S~W 
d 2V 

It will be noted that both the wing and the tail contri- 
butions to dOJd^y are negative; that is, they are in the 
sense to oppose the rotation. 
It is unfortunate that experimental means of measur- 

• rb 

dGJd^y are not more commonly available. As will 

appear later in the report, an accurate knowledge of 
rb * 

dCJd^y is necessary for reasonable accuracy in esti- 
mating stability characteristics. The sparse experi- 
mental evidence concerning the value of this factor (see 
references 7 and 8) indicates that there are, in some 
cases, large interference effects. For one model tested 

on a whirling arm the value of dCJd^for the fuselage 

and tail surfaces combined was only one-third the value 
for the tail surfaces alone. It is quite evident from such 
data that computed values can be considered at best 
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only as rough approximations. With this limitation in 
mind it appears that the most suitable formula for 



dCJd, 



rb 



2 V 



is 



, rb 



3 



l 2 'S t dC L( 



b 2 S dp 



(8) 



and that there is no justification for refinements in the 
formula. 

MASS FACTORS 

Relative density of airplane and air. — The relative 
density of the airplane to the air is usually expressed as 
fi=mjpSb, From this definition p is x/4 times the ratio 
of the mass of the airplane to the mass of air affected 
by a monoplane wing (on the basis of accepted wing 
theory) in traveling a distance equal to the mean chord. 
It thus appears that p is intimately tied up with the 
performance characteristics of the airplane. 

For standard conditions M =-^2^. Yor other con- 



ditions jtt may be expressed as 

lZA(WjS) Po 
b p 



(9) 



where p 0 is the standard mass density (0.002378 slug per 
cubic foot) and p is the actual mass density. It appears 
that jn increases with wing loading and altitude and 
decreases with span. The numerical value of p ranges 
from 2 for large transports to 10 for pursuit airplanes 
under standard conditions. It appears that large air- 
planes are dynamically similar to very lightly loaded 
small airplanes, a transport with a span of 120 feet and 
wing loading of 25 corresponding to an airplane of 
30-foot span with a wing loading of 6.25. 

Ratio of wing span to radius of gyration about X 
axis. — The ratio of the wing span to the radius of 
gyration about the X axis, b/k x , has been determined 
for 15 airplanes (reference 9) and has been found to 
range from 6.7 to 9.3, with 8.0 as an average value. 
This ratio can be estimated with sufficient accuracy for 
stability calculations. For preliminary estimates the 
average value of 8.0 is satisfactory for conventional 
types because, as will appear later, stability character- 
istics are not critically dependent upon the mass dis- 
tribution. 

Ratio of win? span to radius of gyration about Z 
axis. — The value of the ratio of the wing span to the 
radius of gyration about the Z axis, b/k z , has been found 
to vary from 5.1 to 6.4, with 5.7 as an average value. 
As in the case of bjk x , the average value is satisfactory' 
for most estimates of stability. The value of b/k z can 
be estimated with sufficient accuracy for all stability 
calculations from a weight analysis of the airplane. 



STABILITY DERIVATIVES 

In practice it has been found convenient to combine 
the aerodynamic and mass factors that govern lateral- 
stability characteristics into stability derivatives. These 
derivatives take the following forms, one for each of the 
aerodynamic factors: 

1 dC r 

2 dp 

dp 



V*=o 



(£) 



~dp 
2 dCj_ 

d 2V 
dC n 



Change k to k 



_1 (b\ 2 dCn 



'2V 
2 dCi 

jfb 

d 2~V 



d 2V 

Physically these derivatives are, respectively, propor- 
tional to the linear or angular acceleration arising from 
a unit angle of' sideslip, a unit rolling velocity as ex- 
pressed by pb/2V, or a unit yawing velocity as expressed 
by rb/2V. 

The stability derivatives include all the important 
factors governing stability characteristics except p.. 
Since p. occurs only in combination with n g and l v and, 
conversely, since these derivatives occur only in com- 
bination with p, the lateral-stability characteristics can 
be completely expressed in terms of the seven non- 
dimensional quantities: y v , pl v , pn v , l P , n p , l T , and n T , 

For preliminary estimates it will generally be suffi- 
ciently accurate to use the following values for the 
stability derivatives: 

2/,= -0.14 



Z„=32 
w„=16' 



dp 

dCn 

dp 



(10) 



l v =-QA 
n p ~ — 0.5 G L 

U '- 8 ~rb 
d 2V 

A rather small value of y v has been chosen in order to 
be conservative. Stability characteristics calculated 
with this small value of y 0 can be readily corrected to 
correspond to a different ?/„, a fact which will be sub- 
sequently shown. The derivatives l P , n P} and l T may 
differ considerably from the foregoing values, particu- 
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^«»-y)-M»A>0 (12) 

Failure to meet the first of these conditions results in 
spiral divergence," a form of divergence in which the 



larly at angles of attack above that at which the lift- 
curve slope begins to decrease. Fortunately the 
stability characteristics are not greatly affected by 
moderate variation in these particular factors. If 

possible, values of dGJdfi, dCJdfi, and dCJdjpp should 
be obtained by actual measurement. There is strong 
reason for believing that, unless these factors are 
accurately measured, a false impression of the accuracy 
of the estimated stability characteristics may be ob- 
tained by refinements in estimating the other factors. 

FORMULAS FOR ESTIMATING STABILITY 
CHARACTERISTICS 

Stability characteristics about which information is 

fwu precedin S P° rtion of this paper has dealt 

with the various aerodynamic and mass factors that 
govern stability characteristics. In the following para- 
graphs these factors will be grouped in relationships 
which show the effects of the individual factors upon 
the stability characteristics and from which these 
characteristics can be quantitatively determined 

Instability can manifest itself either as a continuously 
increasing divergence from the steady-flight condition 
or as an oscillation of continuously increasing amplitude 
about the steady-flight condition. On a logical basis 
it appears that the questions answered by an estimation 
of stability characteristics should be: (1) Will there exist 

% fwn 11 ?? t0 d ^ er ? 6 fr ° m the Stea dy-Aight condition? 

(2) Will the oscillations started by a disturbance or by 
the use of the controls damp out and, if so, how quickly? 

(3) What will be the period of the lateral oscillations? 
Approximate relationships to answer these questions 
have been developed (see appendixes I and II) and are 
presented in the following pages. 

Basis for formulas.-The following formulas are 
based on the classical theory of small oscillations first 
applied to airplane dynamics by Bryan and developed 
and expanded by Bairstow, Wilson, Glauert, and others 
(references 10 to 13). A brief derivation of the formulas 
is given m appendix I. The formulas presented repre- 
sent a first approximation to a semigraphical method of 
accurately solving the stability biquadratic given by 
the classical theory. This semigraphical method and 
the approximation to it are explained in appendix II 

Formulas for predicting a divergence.— Divergence is 
not possible in the normal-flight range (to which this 
report is confined ; if 

td v n r >nn v l T (id 
and v ' 



airplane tends to go into a spiral dive. Failure to meet 
the second condition results in "directional divergence," 
in which the airplane tends to yaw away from the 
direction of steady flight. 

For purposes of approximate estimation using the 
values for the derivatives given in equation (10) 
equations (11) and (12) become ' 

HdCS®(dCJd^>C L (dC n /d0) (13) 

and 

. -Ciiddidft+zzidCnidftyo (H) 

If the contributions of the wings, the fuselage, and 
interference effects upon dCJd— and upon dCJdfi are 
neglected, equation (13) further simplifies to 



(15) 



This latter equation is, however, an oversimplification 
for any but the most approximate analyses. 
. Formulas for estimating the damping of an oscilla- 
tion.— The number of seconds required for an oscillation 
to damp to one-half its original amplitude is 



(16) 

where f ' is the damping coefficient. The time to damp 
to any other proportion of the original amplitude is 
given by 

where n is the desired proportion, such as % or %. To a 
fairly close approximation (±15 percent) 



s 2 



j Jpn T —l T n v . nn„(— n r ) 
V ^{-n r -l p ) + l-n T -l p y 



+■ 



C L l T n t 



+< 



(~n r ~l p y 



C L l v n T 



21 



-2n,L 



(18) 



In equation (18) the terms in the first pair of brackets 
are those which make f more negative, i. e., decrease 
the time required for the oscillation to damp; the terms 
m the second pair of brackets are those which make 
T less negative. 

If the values from equations (10) are used, equation 
(18) becomes 
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f' = _ 0.07- 



6.4+8 



) 

ia^ ( dC '\ *a dC n /dC n \ 



<-s) 



dC n \ 
jrb 
d 2Vj 



dCn 



(19) 



Since 8 



\ dfi)+cAdfi) 
rE I * s sma ^ com P are( i w ith 6.4 a further 

d h) 

simplification is obtained by letting 8 /—^|-\, where 

\ d 2Vj 

it appears in the denominator, equal 0.84. On this basis 
r = _ 0 .07-3.5^-0.14 <V-1.2 4J 

,/ d<7,\/ „dG. 



+ 0.3 m- 



\ d$rc L d& 



(20) 



Formula (20) will lead to fairly large errors if the air- 
plane departs- very far from the average. The error is 
roughly on a percentage basis so that, for small values 
of damping approaching an undesirable condition, the 
actual error is small. 

Formulas for estimating the period of an oscillation. — 
The period, in seconds, of the lateral oscillation is 



2.83 



(21) 



where \f/ f is the period coefficient. To a fairly close 
approximation \f/' is given by 



Vnh^iip— — nn v l p 
i-n-l v ) 



(22) 



substituting the values* for the derivatives given in 
equations (10) and letting 8^— dO n /d^^=0.84 gives 



P= 



0.146 



_dCA_>^2 dC_ n 

d$rc h d$ 



(23) 



which is correct for the conventional airplane within 
±20 percent. 

CHARTS FOR ESTIMATING STABILITY 
CHARACTERISTICS 

Explanation of charts. — A series of 22 charts for use 
in rapid estimation of stability characteristics are given 
in figures 4 to 25. In these charts the damping and 
the period of lateral oscillations are given by curves of 
Ti-y/W]S^ constant and of P!-y/WjS= constant plotted 
with —{jdCildp as abscissas and tidCJdfi as ordinates. 
The limits to the region within which both spiral and 
static directional stability exist are indicated by straight 
lines representing zero spiral stability and zero direc- 
tional stability, respectively. The rates of convergence 
or divergence are not given. 

vb 

The charts cover values of 6\(0.2 to 2.0) and dC n (d^y 

(—0.030 to —0.252) likely to occur in practice with 
conventional airplanes. 

Each chart covers values of —\LdCi\d$ from 0 to 0.5 
and of ndCJdp from —0.05 to 0.3. These ranges are 
sufficiently large for most conventional airplanes. 
Some extrapolation is permissible in particular cases 
without much loss of accuracy other than that due to 
the fundamental weakness of increasing inaccuracy as 
the damping becomes large. 

These charts are based on equations (13), (14), (19), 
and (23) and are therefore approximations to the same 
extent as the equations. They are intended principally 
for use in rapid estimates in design and show fairly 

vb 

accurately the relative effects of changes in C L , dCJd^yy) 

ndCi/dfi, and ndCJdp. Being based on average values 

i)b 7)b vb 

of dCy/dp, dCi/d~y) dC n /d~y> dCt/d^y* b/k x , and b/k z , 

they cannot be used to determine the effect of changes 
in these factors. The charts should not be used where 
very accurate values are desired. On the other hand, 
there is little justification for using a more accurate 
method unless measured values of the various aerody- 
namic factors are available. If dCyfd^ is known to 
be much larger than —0.28, as, for example, in the case 
of an airplane with a split flap at a high angle of attack, 
correction for the damping can be made by the pro- 
cedure given in the following section. 

Method of using charts. — In order to use the charts 
the following data are needed: 
W/S, wing loading. 
6, wing span. 
C Lf lift coefficient. 
dCild&, rate of change of rolling-moment coefficient 

with sideslip, per radian. 
dCnjdfi, rate of change of yawing-moment coefficient 

with sideslip, per radian. 
An ulk ra ^ e °f °hange of yawing-moment coefficient 
n ' 2V' with rate of yawing, per unit of rb/2V. 
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Figure 9. 

Figures 4 to 9- Lateral-stability charts. 
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Figure 14. Figure 15. 



Figures 10 to 15.— Lateral-stability charts. 
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Figure 17. 





Figures 16 to 21. -Lateral-stability charts. 
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d/3 

l*7S=ieib./sq. ft. 
6=42 ft. CY=0.2 

dCnld—^ 0.060 

Figure 26.— -Special case of lateral-stability chart. 
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From these data, values of —ydd/dp and fidCJdp 
can be readily determined since M = lM^fZg) (standard 
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conditions). 

In general, the value of C L to represent a particular 
range of flight conditions can be chosen as 0.2, 0.8, 1.4, 
or 2.0. It will be necessary in most cases, however, to 

interpolate between two charts for the value of dCJd^y- 
Any point given by -ndCt/dp, pdCJdp represents a 
value of Pf-y/W/S and a value of T/^W/S. The 
period and the time to damp to one-half amplitude are 
readily obtained by multiplying these values by -y/W/S. 
The location of the point -pdCt/dp, fidCJdp also 
indicates whether there will be a tendency to diverge. 

The charts are computed for standard conditions. 
They can be easily applied to a study of stability at alti- 
tude by substituting the value (0.00238/p) W/S for 
W/S wherever W/S occurs in the computations. 

Correction to a different value of dC Y /dp may be 
readily made as follows: Compute £' ioirdO Y /dp== —0.28, 



i. e., y t =— 0.14, from the relation f 



,_ -0.313VW/S 



Add to this value of f ' the quantity ~(0.2S+dC r /dp) to 

obtain the co rrect ed value of f . Calculate the corrected 
value of T/tJW/S using the corrected value of f'. 

In cases where a large number of estimates are to be 
made for a given pair of values of W/S and b, it will 
sometimes be convenient to convert the charts to read 
directly in terms of -dQ/dp, dCJdp, P, and T. This 
conversion can readily be accomplished without re- 
drawing the chart by changing the constants. Figure 26 
represents figure 5 converted to read directly in the 
desired quantities for an airplane having W/S=1Q 
pounds per square foot and 6=42 feet. 

Example of use of charts.— It is assumed that the 
lateral-stability characteristics throughout the normal- 
flight range are desired for a 5,000-pound airplane having 
a wing loading of 16 pounds per square foot and a span 
of 42 feet. Values of dCJdf and dCJd* are available 



from wind-tunnel tests. 



rh 

Values of dCJd— must be 

estimated. The airplane is a modern type with a fairly 
high top speed and is equipped with split flaps. Flaps 
were considered to be down at <7 L =2.0 but up at 
O L =0.2, 0.8, and 1.4. 

The stability characteristics will be estimated for 
each of the C h values of 0.2, 0.8, 1.4, and 2.0. Values of 
—tidC if dp and ndCJdp are determined at each value of 
C L from the relationships 



and 



■^dC l /d^^^X57.3X(dC l m) 



f >dCJdp= 1 -^p*X57.3X(-dC n /d*) 



Values of dCJd—^&Te determined from the relationship 

where C Dw is taken from wind-tunnel tests of a similar 
wing, I /b and SJS are dimensional characteristics of 
the airplane, and dC ht /dP is estimated using the 
relationship 

dC Lt _ 5.5 



dp 



1 



h 2 /S t 



where b t is the height of the vertical tail surface. The 
values of C Dw , dCJd~, -vdC t /dp, and ixdGJdp are as 
follows: 



Ct 




4<W& 


-ndCldfi 


ndCn/dff 


0.2 
.8 
1.4 
2.0 


0.008 
.025 
.070 
.400 


-0.051 
-.057 
-.072 
-.182 


0.20 
.25 
.45 
.55 


0.180 
.180 
. 165 
.130 



From the various charts, values of T and P are de- 
termined, interpolations and extrapolations being made 
where necessary. The values of the stability charac- 
teristics at each value of C L follow. 



Ct. 


T 
(sec.) 


P 
(sec.) 


Divergence 


0.2 
.8 
1.4 
2.0 


2.1 
3.1 
3. 1 
1 2.2 


3.0 
5.5 
5.8 
6.4 


None. 
Spiral. 
Do. 
None. 



va'lfeT r^lSsecondr 6386 dCrW dm k ° the high drag gives the correeted 

EFFECT OF THE GOVERNING FACTORS ON THE 
STABILITY CHARACTERISTICS 

AERODYNAMIC FACTORS 

lateral force due to sideslip.— The lateral force due 
to sideslip is small, in general, but beneficial in its 
effect upon stability characteristics. As appears 
in equation (18), dC r /dp adds directly to the damp- 
ing coefficient, A£' = y 4 AdC y /dp. For the value of 
dC Y /dp= -0.28, Af'= -0.07, which is sufficient to 
damp the lateral oscillation to one-half amplitude in 
8 seconds for an airplane with a wing loading of 16 
flying at 176 miles per hour. The effects of dC r /dp 
on the period and on the tendency to diverge are 
negligible. 

Rolling moment due to sideslip.— The rate of change 
of rolling-moment coefficient with sideslip plays a great 
part in determining the stability characteristics, as is 
apparent from a glance at the charts of figures 4 to 
25. It is necessary for stability that dC t /dp be negative ; 
the term —dG % /dp will be used, as in the charts, for 
simplicity in discussion. 
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Increasing —dCJdft. increases the range of values of 
dCn/dp within which there is no divergence, there being 
less likelihood of either spiral divergence or directional 
divergence as —dCt/dp is increased. Increasing —dCi/d$ 
increases the time for an oscillation to damp and short- 
ens the period. These effects are sufficiently small to 
be of no practical importance at high speeds but are 
appreciable at low speeds. 

Yawing moment due to sideslip. — From the con- 
siderations of tendency toward divergence the value of 
dCJdp should be small and positive. Too large a 
positive value of dCJdp results in spiral divergence. 
Too large a negative value results in directional 
divergence. The range of permissible values of dCJdp 
is quite narrow for small values of —dCi/dfi and 

— dCnjd^- 

Increasing dCJdfi increases the damping and shortens 
the period of the lateral oscillations. The effect upon 
the period is very pronounced, particularly at small 
values of the lift coefficient corresponding to cruising 
and high speeds, as is especially apparent in equation (23) 
where 

Olib 




{-dCtldft+ykcjdP 

It appears that for 0^=0.2 the effect of dCJd& upon 

the period has 16 times the effect of —dCi/d(3. The 

effect upon T is less pronounced. It is of interest to 

note that the theory indicates stability with dCJdfi 

• vb 
zero or slightly negative if dC Y /d$ and dCn/d^y&ve 

moderate or large. 

Rolling moment due to rolling. — Differences in the 

value of dCijd^yoi the order of those likely to exist 

between conventional airplanes in the normal-flight 
range have but slight effect upon the tendency toward 
divergence or the oscillatory characteristics. This fact 
tends to justify the use of an average value for this 
factor in equations (20) and (23) and in the charts. 
The small effects occurring are such that increasing 

dCi/d^- decreases the time required to damp, in 

general, and increases the period. 

Near the stall dO t /d^ changes sign and tends to 

result in violent instability. This report does not deal 
with stability near the stall, which is amply discussed 
in references 14, 15, and 16. 

Yawing moment due to rolling. — As is the case for 

dCi/d^y) differences in dCJd^r likely to exist in 

practice have comparatively slight effect upon the 
stability characteristics below the stall. Increasing 

dCJd^y may either increase or decrease T, depending 



upon the magnitudes of other quantities, and increases 
P slightly. Here again the selection of an average 
value for this factor seems justified. Near the stall 

dC n {d~y changes sign and becomes an important 

factor in producing instability. 

Rolling moment due to yawing. — The rolling moment 
due to yawing is chiefly of importance in connection 
with the likelihood of spiral divergence. Increasing 

dCi/d-^y decreases the range of values of dCJdfi for 

which spiral convergence exists for a given set of 

vb vb 
values of —dCi/dfi and dCJd^y- Increasing dCi/d^y 

generally decreases T but has no noticeable effect upon 
P or the likelihood of directional divergence. 

vb 

Yawing moment due to yawing. — Increasing dCJd^y 

increases the permissible range of values of dCJdfi for 
spiral convergence and decreases the time required for 
the oscillation to damp to one-half amplitude. It is 
apparent from equations (20) and (13) and from the 

vb 

charts that an accurate knowledge of dCJd^y is 

essential to accurate calculations of T and of the 
limiting values of dCJdfi within which spiral conver- 

vb 

gence exists. On the other hand, dCJd^y has only a 

very slight effect upon the directional convergence 
or upon the period of the oscillations. 

MASS FACTORS 

Relative density of airplane to air.— The relative 
density n has no effect upon the likelihood of either 
spiral or directional convergence. Its effect upon the 
period and damping of the lateral oscillation can best 
be understood by considering the separate effects of the 
factors which determine p, namely W/S, b, and p. 
Since a decrease in p has precisely the same effect as 
an increase in W/S, the effects of altitude are the same 
as the effects of increasing the wing loading and will 
therefore not be discussed separately. 

The effect of wing loading upon the time required to 
damp the oscillation to one-half amplitude can best be 
deduced from equations (16) and (20). From equa- 
tion (20) it appears that if, for the case at hand 

dO, / dC n 

dp 



1.2 



(-ft) 



is greater than 0.3 (—dCJdp), then 



increasing W/S (since p is proportional to W/S) will 
make f greater in the negative sense. This will be the 
case only for very small values of —dCi/dfi. In gen- 
eral, therefore, increasing W/S will increase T both by 
decreasing f' and by increasing the numerator in the 
relationship 

T -0.313 -yj(W/S)C L 
2 ~ f 
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On the other hand, wing loading has no appreciable 
effect upon the period, at a given C L , as is apparent from 
equation (23). 

From the charts it appears that, since increasing b 
decreases ndCJdfi and —ndCJdp, increasing the span 
will decrease the time required to damp the oscillations 
over most of the range of values of the parameters. As 
pointed out in the preceding paragraph, the effect de- 
pends upon the relative magnitudes of dC n /d0, dCJd^y, 

and —dCi/dp. Only in the case of very small values of 
—dCi/dfi, will increasing the span increase T. For prac- 
tical purposes the period of the lateral oscillation is 
proportional to the square root of the span, as is shown 
by equation (23). 

Ratio of wing span to radius of gyration about X 
axis. — In the discussion of the effects of changing 
b it was assumed that the ratios b/k x and b/k z were kept 
constant. The effects of changing these ratios can be 
most readily explained on the basis of keeping b 
constant. 

The value of k x has no effect upon either spiral or 
directional convergence. Although not readily appar- 
ent in equations (18) and (22), increasing k x results in 
small increases in T and P. There is, however,, no 
justification for extensive labor to determine k x accu- 
rately in the absence of accurate data on all the aero- 
dynamic factors. 

Katio of wing span to radius of gyration about Z 
axis. — The effects of increasing k z are similar to the 
effects of increasing k x . It has a slight but unimportant 
effect sUpon directional convergence. Its effect upon 
the period is greater than the effect of increasing k X) but 
not great enough to be of practical importance in most 
cases. 

GENERAL COMMENTS 

The present state of knowledge does not justify 
positive assertions as to the desirability of any given 
set of stability characteristics. Very little has been 
done to determine quantitatively the stability charac- 
teristics that result in the most satisfactory riding and 
handling characteristics. Such research (reference 17) 
has given more or less negative results, at least with 
respect to the period and the damping of oscillations. 
It is definitely known, however, that very great instabil- 
ity, such as that at the stall, and very great stability 
are both undesirable. There is strong reason to believe 
that any tendency to diverge i« undesirable but that, 
if such a tendency is of small magnitude, it will not 
seriously inconvenience the pilot. 

When the foregoing facts are taken into considera- 
tion, it seems desirable that for airplanes designed for 
most purposes, excepting machines intended as pur- 
suits, fighters, or for acrobatics, there should be no 
tendency to diverge, oscillations should be moderately 
to heavily damped, and the period of the oscillations 
should be as long as practicable. It is believed that 



such characteristics will require a minimum of effort 
from the pilot and will result in a maximum of passenger 
comfort. 

Reference to the charts of figures 4 to 25 reveals that 
these characteristics can be attained only by making 

dCJdp small while keeping dCJd^. and dC Y fdp large. 

Some additional advantage is gained by keeping 
—dCi/dp small, particularly at high angles of attack. 
Probably the best value of —dC t [dp is from small to 
moderate, the moderate values giving more pronounced 
spiral and directional convergence. The best method 
of keeping dC n /dp small while retaining large values of 
vb 

dCJd^y and dC Y \d$ appears to be the use of a fuselage 

giving an unstable yawing moment of rather large 
magnitude. The unstabilizing effect of the fuselage 
depends on its length, breadth, the distance of the 
center of gravity from the nose, and on the shape. 
The shape of the fuselage, and possibly interference 
effects, play an important part, which can be deter- 
mined accurately only by wind-tunnel tests. A large 

rb 

value of I tends to make dCJd^y large and a large, 

deep fuselage tends to give a large value of dC Y /dp. 
The dihedral of the wings can be adjusted to bring 
dCi/dfi to the desired value but here again, with present 
knowledge, it is necessary to make wind-tunnel tests. 

The formulas and charts indicate that spiral diver- 
gence is likely to occur at slow speeds and that direc- 
tional divergence may occur at high speeds. Increas- 
ing oscillations do not appear to be possible within the 
normal range of aerodynamic governing factors. In 
some Cases, however, the damping is so slight that the 
oscillations might be unpleasantly apparent. 

Although control is outside the intended scope of 
this paper, it should be pointed out that appearance of 
instability may, under certain circumstances, be brought 
about by the influence of the controls. The two most 
common instances are that of directional divergence 
arising out of an attempt to hold the wings level with 
conventional ailerons, the rudder being held neutral; 
and that of increasing or poorly damped oscillations 
arising out of operation of the rudder in improper 
phase relationship to the change in attitude of the 
airplane. The directional divergence is caused by the 
adverse yaw of the ailerons and can be avoided by re- 
ducing the adverse yaw, by increasing dC„/d@, or by 
holding the ailerons neutral and allowing the airplane 
to roll. The increasing oscillations are most likely to 
occur when the natural period of the airplane is short 
and when the rudder is operated in such a manner as to 
prevent yawing. They can be avoided by holding the 
rudder neutral or by operating it in such a manner as 
to produce sideslip opposing the roll, i. e., by trying to 
hold the wings level rather than by trying to prevent 
yawing. 
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SUGGESTIONS FOR FUTURE STUDY 

A systematic correlation of stability characteristics 
with riding and handling qualities is needed. It is 
possible that these qualities are more directly related 
to certain of the governing factors than to the tendency 
to diverge or to the characteristics of the oscillations; 
investigations should be conducted with this possi- 
bility in mind. - 

There is need for adequate comparison between com- 
puted values and measured values of stability char- 
acteristics as a check upon the accuracy and validity 
of the mathematical treatment. 

At present some of the aerodynamic governing factors 
cannot be estimated with assurance. A great deal of 
systematic study will be necessary to provide sufficient 
data for the formulation of satisfactory empirical con- 
stants to be used in estimating these factors. 



In this report no consideration has been given to the 
effect of power on the lateral stability. This problem 
should be the subject of a study sufficiently thorough to 
reveal the effects of power on the stability derivatives 
and upon the mathematical treatment necessary to 
estimate the stability characteristics. 

More satisfactory means of measuring the separate 
aerodynamic factors and also the final stability char- 
acteristics of models are necessary for rapid progress. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., November 17, 1936. 
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DERIVATION OF FORMULAS 

The theory of small oscillations. — The theory of 
small oscillations was first applied by Bryan to the 
dynamics of mechanical flight (reference 10). On the 
assumption that the direction and magnitude of changes 
in the aerodynamic characteristics due to changes in 
motion from the steady-flight condition are known? 
equations of motion in unsteady flight are written for 
the case of small deviations from the steady condition, 
one equation for each of the degrees of freedom of the 
motion. Simultaneous solution of the equations gives 
values that describe the motion of the airplane after a 
disturbance. 

Assumptions in the application of the theory.— In 
the application of the theory of small oscillations to 
quantitative estimations of stability characteristics, a 
number of assumptions are necessary in order that the 
mathematics may not be too involved and the compu- 
tations too extensive for practical applications. The 
primary assumptions are as follows: 

(a) The combined aerodynamic effect of two or more 
components of motion is assumed equal to the alge- 
braic sum of the separate effects of the individual 
components. 

(b) The changes hi aerodynamic forces and moments 
due to a deviation are assumed proportional to the 

deviation, i. e., the slopes dCt/dp, dCtfd^yt etc., are 

assumed to be constants. 

(c) The lateral motion involving p f q, and r is as- 
sumed to be independent of the longitudinal motion, i. e., 
the machine is assumed to be symmetrical. 

(d) Secondary effects such as those involving the 
products of two or more small quantities are neglected. 

(e) The values of the aerodynamic factors are 
assumed to be unaffected by the linear and angular 
accelerations. 

Equations of lateral motion. — The equations of 
lateral motion will be written for the axes shown in 
figure 1 using the symbols and notation given in 
appendix III and on the report covers. The X axis 
is taken in the direction of the relative wind during 
the steady-flight condition. The axes are assumed 
fixed in the airplane. During steady flight, 

Y=L=N=0 
v=p=r=0 

After a disturbance, 



dY, dY, dY, m . . 
v W + ^+ r W+ W sin * cos y 



+ W^sin ^ sin y=m^-\-mru 



dL . dL , dL , 9 uru 
v Tv+*di+ r ^ mk *ft 



dN , dN , dN 



dr 



dv' F dp 



dr 

p—d4/dt 
r=dif//dt 



(24) 



(25) 



It is assumed in these equations that the principal 
axes of inertia are coincident with the reference axes, 
which is not true in the general case. A number of 
supplementary calculations made as part of the study 
leading up to this report have indicated, however, 
that to neglect the angularity of the principal axes to 
the reference axes will not introduce serious error in 
the normal-flight range and will give slightly conserv- 
ative results. Consequently, the terms including the 
product of inertia were omitted to make the equations 
as simple as possible. 

Since dY/dp and dY/dr are small, they are generally 
neglected. For the small deviations considered, u may 
be taken equal to the steady-flight velocity V and the 
sines of the angles of roll and yaw may be replaced by 
the angles themselves. Since in power-off flight the 
lift is equal to W cos y and the lift times the tangent of 
the angle of glide is equal to Wsin y, the first of the fore- 
going equations will be rewritten, 



dY 



dv 



0^+0X(Uft) + *X(Uft)Xtan y = m %+mrV 



dt 



The equations of equilibrium finally become, 



dY dv 



dv -™^+0Oift) + *aift) t&n y—f t mV=0 
d\p dL 



dL . d(j) dL d 2 <f> 
dv dt dp dt' 



dv 

Replacing 



dN djdN ,d*dN_W, 
'-■ + dt dp^ dt. dr dt 2 miCz ~ U 



(26) 



dCy 

f 

(Lift) by£pVW<7* 
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% by 1,™^ etc. 

fb y ^ 6 f»,etc. 
m , 



by 



dv Dy V 



Writing in determinant form and simplifying gives, 
vy °~ T dt T* 



§r (tony) ip—r 



dt 



v id,, 
v nn, 



lpT dt dt 2 



lfT dt 



n » r Tt 
where V'=kw 

' 'Akz) d& 
j^l/bVdGt 



UrT dt 



=0 (27) 



> 



2 dCi 



d 



rb 

W 



Substituting y 0 e>< for 0 O c xl for 0, etc., and simplify- 
ing give 



pi* 



*2* *2~ tan 7— rA 



Z p tX-t 2 X 2 
7l„rX 



/ r rX 
71 t t\-t 2 \ 2 



=0 



(28) 



which can be expressed as 

A(\y+B(yy+c(\y+D(\y+E(\')=o (29) 

• where, 
X'=rX 

B=—yt>—n r —l p 

C= l p n r — Z r 7i P + y v (l p + n r ) + nn v 

D=yi(kb-lfr)+rfj(n,-^-HnJ^ tan 7 +Z p ) 



•i r n t )+ii 



C, 



2 ^^r-^^T/fy tan ~i(l v n v —l 9 n v ) 

The solution X'=0 is readily apparent in equation 
(29) and results from the fact that the airplane has no 
inherent tendency to return to any fixed compass course, 
being, as the solution shows, neutrally stable in that 
respect. This solution, X'=0, is generally neglected 
and the lateral-stability characteristics are considered 
as those given by the biquadratic 

' AQi'y+B(\y+CW+D(\') + E=0 (30) 
The deviation of each one of the components of the 
lateral motion varies with time according to the rela- 
tion, 

v, p, or r^de^+G^+C^' + G^** ' 
where X x , X 2 , X 3 , and X 4 are the four roots of the bi- 
quadratic. In lateral motion the constants B, G, D, 
and E are generally such that there are one pair of real 
roots and one pair of conjugate complex roots indicating 
motion of the type, 

v, p, orT=(7 5 ^(cos W-Cd+CW+CnJ* 
where f and ^ are the real and imaginary parts, respec- 
tively, of the conjugate roots. Tins motion represents 
an oscillation superimposed upon two rates of con- 
vergence (or divergence) . It is evident that for stability 
X 3 , and X 4 must be real and negative so that the 
values of v, p, and r will reduce to zero. In order that 
the real parts of the roots shall all be negative it is 
necessary and sufficient that B, C, D, E, and 
(BCD~D 2 -B 2 D) each be positive. In order to 
determine the rates of convergence and the damping 
and period of the oscillation, it is necessary to solve the 
biquadratic. A convenient semigraphical method of 
solving stability biquadratics was pointed out in 
reference 1 and is described in detail in appendix II. 



APPENDIX II 



SOLUTION OP STABILITY BIQUADRATIC 

Semigraphical method of solving biquadratics.— The 
biquadratic 

X 4 +#X 3 +C7X 2 +Z}\+#=0 
can be expressed as 

(\>+a x \+b x ) (x 2 +a 2 X+6 2 ) = 0 

from which 



(31) 



(32) 



are roots of the general equation. It appears that 
2?=a 1 +a 2 

C=a 1 0 2 + &1 + &2 

D=a,ib 2 +a 2 bi 
E==b x b 2 

Eliminating values of and b 2 , 



and 



_ b x 2 B-b x D 
b x 2 -E 



(33) 

(34) 
(35) 



Note that, if the minus sign in equation (34) is chosen 
for a u the plus sign will correspond to o^. 

Values of a x and b x that wul satisfy these equations 
separately are plotted on charts having values of a 
as abscissas and values of b as ordinates. The inter- 
section of the resulting curves represents values of a x 
and b x that satisfy both equations. There are two 
intersections in the general case, one corresponding to 
a x and b x , the other to a 2 and b 2 . Ordinarily it is more 
convenient to find one of the intersections by plotting 
and to solve for the remaining values by the use of 
equations (33). 

Figure 27 gives a solution of a typical stability bi- 
quadratic and illustrates the use of this semigraphical 
method. For most cases time can be saved in locating 
the intersection by letting a x equal zero in equation (35), 
thus dete rminin g the intersection of the plot of equa- 
tion (35) with the b axis. The resulting value of b x 
when substituted in equation (34) will give an approx- 
imate value for a x . The final values can then be de- 
termined as accurately as desired by locating the point 
of intersection of the curves. This method has been 
applied to several hundred solutions of stability bi- 
quadratics in the course of the study leading up to this 



report and has been found to be very satisfactory, 
particularly so if systematic changes in factors are being 
studied. 



-.2 



i I i r 




J 5.52 -t //5.52\2 ~ _ ~~ 
— -V(— )-5.36 + 6 l + 



0.74 



.6 



1.0 



Figure 27.— A semigraphical method of solving stability biquadratics. Solution of 
biquadratic: 

X«+5.52\3+5.36XM-13.90X+O.74=(X'+O|X+6,)(X>+aj\+6j)=0 
From intersection of curves 

ai-0.46; 6i=2.72 

From equations (33) 



aj=5.52-0.46=5.06; 6j=0.74/2.72=0.027 



From X»+0.46X+2.72 =»0, X= -0.23±i 1.6 
From X»+5.06X+0.027=*0, X=-5.05 
X— 0.01 



Approximate formulas for the dampinj^nbV'tJie* 
period of the oscillation. — As was stated in the preceding 
paragraph, an approximate value of bi can be found by 
substituting a x =0 in equation (35) giving 



i D 

1 F 



(36) 



Substitution of this value of b x in equation (34) gives 
the approximate value for a x of 

This latter equation can be further simplified without 
loss of accuracy by removing the radical and assuming 
a x 2 =0. The resulting equation for a x is 



Since 



ai ~B B* D 



^-Wd)'-*- 

s 2 



(38) 

(39) 
(40) 
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The period in seconds of the lateral oscillation is 



*' = V-(f)' +6 ' 



Since '(jjfy is normally small compared with b u it is 
sufficiently accurate for practical purposes to put 

*'=V£ (41) 
From equations (38) and (36), 

Id 



(43) 



Supplementary study has indicated that further simpli- 
fication can be had with but slight loss in accuracy by 
neglecting values of y v in the expression for B and D 
and neglecting tan y in the expression for D and E 
giving 

B=—n r —l p 

C=l P n r ~lrn P +y,(l P +n r )-j-iin n 
■D=tfj(n 9 -Q)-Mjl p (44) 

E=n~ (lvn r ~l T n^ 

Substituting these values in equations (42) and (43) 
and simplifying where possible gives 




l p n T —l T n P , nn v {~ n t ) 



{—Thr—k) 



{~n r -l p y 
0 L l 



2l t (n } 



f)-2*4 



+1 



tf 9 (n p —~^ 



CJ, v n r 



(~n r -l v ) 2 21 



and 



V (-n-l P ) 



v(n P -~y 



■2nJ, P 



(45) 



(46) 



The time in seconds to damp to one-half amplitude is 
given by 

T_Jog« 0-5 -0.693 
T — p— r 

(smce t'=T£). 

Expressed in more convenient form, this equation 
becomes 



r 



(47) 



or 



p== 2T = 2jrr 



(48) 



Approximate formulas for the convergence char- 
acteristics. — Since ai is small, a 2 is approximately equal 
to B. 

Letting a 2 =B 
and 



7 EB 



gives 



(49) 



as the solution for the pair of the roots of the stability 
biquadratics corresponding to the convergence charac- 
teristics. Since B is always positive and large in the 
normal-flight range, it appears that this equation repre- 
sents two convergences if EB/D is positive and less than 
(B\2 ) 2 , a heavily damped oscillation if EB/D is positive 
and greater than (B/2 ) 2 , a divergence and a convergence 
if EB/D is negative and less than (B/2) 2 , and two di- 
vergences if EB/D is negative and greater than (-B/2) 2 . 
Instability is therefore possible if either E or D be- 
comes negative. For most cases E is small but may 
be either positive or negative and D is positive and 
large. These circumstances give the usual solution 
of (49) as a large negative root approximately equal 
to ~B—n T -\-l p and a small root approximately equal 
to -E/D. 

In the usual case it is desired to know whether or 
not there will be a divergence rather than to know the 
rapidity of the convergence. For such a case it is 
sufficient to know that 



and 



Z>>0 
E>0 



By the use of the relationships of equations (44), these 
conditions are represented by 



I 



>(n P — y)~ M>>° 



and 



hn T y>l r n v 



(50) 



(51) 



These equations neglect the effects of y e and tan y, a 
procedure that is conservative for power-off flight. 



TV 



APPENDIX HI 



SYMBOLS 

X, Y, Z, axes of reference fixed in the airplane having 
the origin at the center of gravity, the X 
axis in the plane of symmetry and along 
the relative wind in steady flight, the T 
axis perpendicular to the plane of sym- 
metry, and the Z axis in the plane of 
symmetry and perpendicular to the X axis. 
X, Y, Z, forces along the respective axes, X being 
positive when directed forward, Y positive 
when directed to the right, and Z positive 
when directed downward. 
L t M, N, moments about the X, Y, and Z axes, respec- 
tively, L being positive when it tends to 
depress the right wing, M positive when it 
tends to depress the tail, and N positive 
when it tends to retard the right wing, 
components of linear velocity of the airplane 
along the X, Y, and Z axes, respectively, 
having the same positive directions as the 
X, Y, and Z forces, 
resultant velocity. 

components of angular velocity about the X, 
Y, and Z axes, respectively, having the 
same positive directions as L, M, and N. 
4>, 0, rp, components of angular displacement from a 
given attitude about the X, Y, and Z 
axes, respectively, 
a, angle between the relative wind on a plane 
parallel to the plane of symmetry and the 
wing chord, positive when corresponding 
to positive rotation 6 of the airplane 
relative to the wind. 
jS, angle between the relative wind and a plane 
parallel to the plane of symmetry, equal to 

sin _1 p.; angle of sideslip in radians. 

7, angle of flight path to horizontal, positive in 

a climb. 
T, dihedral angle, degrees. 
Y 

C r =-o, coefficient of lateral force. 
qS' 

C,=-or, coefficient of rolling moment. 



V, 

P, 2, r, 



^ n= qM' coem °* ent °* yawm S moment. 



ft, 



coefficient of drag for the wing alone, 
coefficient of force on projected side area, 
coefficient of force on vertical-tail area. 



S, wing area. 

S„ projected side area of fuselage. 
S t , vertical-tail area. 
Sp, projected side area. 

d, maximum depth of fuselage. 

y, spanwise distance from plane of symmetry. 
x u distance from fuselage nose to center of 
gravity. 

I, distance from center of gravity to rudder 
hinge. 

li, over-all length. 

Z 2 , over-all length of fuselage. 
Po, mass density of air under standard condi- 
tions. 

p, mass density of air under condition of flight. 

t, subscript denoting vertical tail surfaces. 
b t , height of vertical tail. 
rj h tail efficiency. 

z p , the Z coordinate of the center of pressure of 
projected side area. 
K#, empirical factor for estimating dC n /dp for 
fuselage. 

m 



m=-ot, relative density factor. 



P Sb 



For standard atmosphere, ju= 



iz.UW/S) 



T=-^ r7 =n~, time conversion factor. 
poK V 

IdC 

y»=2 ~dB ' non(umens i ona l derivative of lateral force due 
to sideslip. 



l dCi 



> nondimensional derivative of rolling 



moment due to sideslip. 



' 2\kJ 



1/6 \ 2 dC 

n9== 2\F~) ~di' noncumensionai Privative of yawing 
^ ^ moment due to sideslip. 

^iCfT^T' nondimensional derivative of rolling 
** d$rf moment due to rolling. 



2V 



Wp =I^Y^^, nondimensional derivative of yawing 

d 2V moment due t0 rollm g- 
; r =l/AY^., nondimensional derivative of rolling 
d W moment due t0 vawm S- 

r =l^Y^, nondimensional derivative of yawing 
* K 2 J - } - moment due to yawing. 



d 2V 

B, C, D, E, coeflicients of stability biquadratic. 
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^=~=tdti^=—±i—f root of stability equation. 

T, time for oscillation to decrease to one- 
half amplitude, seconds. 
-0.693 K , 

■AT 



6. 



= zM^mm (Standard atmos- 
* phere). 
P, period of oscillation, seconds. 

-A' 

j^jmmou (Standard - atmos . 

^ phere). 
a, 6, coefficients of stability quadratic. 
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